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Abstract: This paper presents an optimum design of a lightweight vehicle levitation 
electromagnet, which also provides a passive guide force in a magnetic levitation system 
for contactless delivery applications. The split alignment of C-shaped electromagnets about 
C-shaped rails has a bad effect on the lateral deviation force, therefore, no-split positioning 
of electromagnets is better for lateral performance. This is verified by simulations and 
experiments. This paper presents a statistically optimized design with a high number of the 
design variables to reduce the weight of the electromagnet under the constraint of normal 
force using response surface methodology (RSM) and the kriging interpolation method. 2D 
and 3D magnetostatic analysis of the electromagnet are performed using ANSYS. The 
most effective design variables are extracted by a Pareto chart. The most desirable set is 
determined and the influence of each design variable on the objective function can be 
obtained. The generalized reduced gradient (GRG) algorithm is adopted in the kriging 
model. This paper's procedure is validated by a comparison between experimental and 
calculation results, which shows that the predicted performance of the electromagnet 
designed by RSM is in good agreement with the simulation results. 

Keywords: optimum design; response surface methodology (RSM); levitation control; 
kriging interpolation method; design of experiments (DOE); electromagnet 
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1. Introduction 

Electromagnetically levitated and guided systems with linear motor propulsion are commonly used in 
the field of people transport vehicles, tool machines and conveyor systems because of their silent and 
non-contact motion [1-3]. Passive guidance controls are normally used in real implementations of 
Maglev systems with linear induction motor propulsion to reduce construction costs and complexity. The 
Japanese HSST and the English BAMS (Birmingham Airport Maglev System) are well-known cost- 
effective transport systems based on linear induction motor propulsion and magnetic levitation. In both 
systems, the guidance force needed to keep the vehicles on the track is obtained by the levitation 
electromagnets with the help of particular shapes of the rail and clever placement of the electromagnets 
about it. In Maglev transport systems where levitation magnets are on moving parts, the design of 
levitation electromagnet is the most important factor. It can account for a considerable part (10%) of the 
whole vehicle weight, which can largely affect weight and stability of the magnetic levitation system. 
Optimum plans for the weight reduction of the levitation electromagnet are consequently 
desirable [4-7]. In this paper, an optimal design procedure was required in order to develop a magnetic 
levitation system for a contactless delivery application shown in Figure 1 . As a first step, the design 
goal was to reduce the weight of the electromagnet of the magnetic levitation system with the 
constraint of normal force, considering the initial model. At a second step, the most effective design 
variables and their levels should be determined and be arranged in an orthogonal array table. 

Figure 1. Magnetic levitation system prototype for contactless delivery application: (a) 3D 
modeling (b) prototype. 



Response Surface Methodology (RSM) is generally used with two or three design variables, 
however we have seven design variables, so a mixed orthogonal array table was utilized. The reason 
for this is that a mixed orthogonal array table is an efficient way to study the effect of several design 
variables simultaneously with a number small of experiments and to plan matrix simulation trials. For 
each design variable combination the response value is determined by the 2D and 3D Finite Element 
Method (FEM). In this paper we use the reduced gradient algorithm, which can lead to the selection of 
the most desired set of variables. Figure 2 shows optimization procedure flow chart of RSM and 
kriging interpolation method. Figure 3 shows the application prototype model for applying the 
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contactless delivery application. The response is determined and we evaluate the influence of each 
design variable on the objective function. Based on this method the weight of the optimized 
electromagnet could be reduced by 11.412% and normal force improved by 7.754%, compared to the 
initially designed electromagnet. 



Figure 2. Optimization flow chart. 
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Figure 3. Application prototype. 




2. Passive Guidance Control and Optimization of Electromagnet 



Let's assume that the guidance forces are generated by closed-loop control of levitation 
electromagnets. The lateral response of the electromagnet due to the airgap control is a force 
increasing with the lateral offset X. This force is almost a linear function of the lateral offset of the 
electromagnets. Then we can imagine that the electromagnet will move laterally as a mass bound by a 
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spring. As we know, a spring working on a mass is a mechanical resonant system, whose resonance 
frequency is given by: 



f resonance r\ 




(1) 



where the M is the sum of the masses of the electromagnet and of part of the carried vehicle. We can 
easily imagine that any external action will cause a non-damped oscillating response of the lateral 
position. The value of k gu id e can be designed to keep the disturbance from causing an excessive 
displacement. An example of the design problem is shown in Figure 4. 

Figure 4. Lateral position model for electromagnets under constant levitation control: 
(a) passive guidance model; (b) simulation result for a 100 N disturbance. 
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A total mass of 100 kg is assumed to be levitated with a constant gap. The sum of lateral position 
stiffness of the levitation magnets is set at 10,000 N/m, 50,000 N/m and 100,000 N/m, respectively. 
The lateral disturbance forces are assumed to be a constant value of 100 N for about 1 second. The 
simulation result shows that if the lateral position stiffness remains higher than 50,000 N/m, the lateral 
position deviation is smaller than 5 mm under constant 100 N disturbance. 

In the experiment, a total mass of 200 kg including four levitation electromagnets is levitated under 
small deviation of electromagnet placement conditions, as shown in Figure 5. That is, a mass of 50 kg 
is levitated by one electromagnet. In the experiment, levitation magnets are controlled to maintain a 
constant gap length of 5 mm. A position disturbance of 3.2 mm is applied at first and then removed. 
The lateral response of the experiment is shown in Figure 5(b). We can infer that if two electromagnets 
are positioned with small deviation from the original rail position, the guidance forces generated each 
electromagnet can be larger, but the resultant force becomes smaller because of the differential 
actuation scheme, so the position stiffness become smaller, the natural frequency becomes smaller and 
the lateral deviation becomes larger as the deviation of magnets becomes larger. Therefore, no-split 
positioning of electromagnets is better for lateral performance. 
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Figure 5. Lateral position response of magnetically sprung mass under position 
disturbance: (a) experimental model; (b) lateral position response (experiment); (c) lateral 
position stiffness response (simulation); (d) lateral force response (simulation). 
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3. Optimum Theory 

3.1. Response Surface Methodology 

The RSM method can be readily adapted to develop an analytical model for a complex problem. 
With this analytical model, an objective function with constraints can be easily created and evaluated, 
and computation time can be saved. A polynomial approximation model is commonly used for a 
second-order fitted response (u) and can be written as follows: 

u = (3 0 + £ P,x, +i^x iXj +s (2) 

H j=i j=l 

P : regression coefficients, x : design variables 
S : random error, k : number of design variables 

The least squares method is used to estimate unknown coefficients. Matrix notations of the fitted 
coefficients and the fitted response model should be as shown in equation (3): 

$ = (XXy l X'u 9 u = J$ (3) 
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where P is a vector of the unknown coefficients which are usually estimated to minimize the sum of 
the squares of the error term. It should be evaluated at the data points. RSM can be applied in 
connection with FEM and the response actually represents FEM result. 

3.2. Kriging Interpolation Method 

Kriging is a method of interpolation named after a South African mining engineer D. G. Krige, who 
developed the technique while trying to increase accuracy in predicting ore reserves. In the kriging 
model, the global approximation model for a response y(x) is represented as: 

y(x)=j3+v(x) (4) 

where x is the design variable vector, /? is a constant, and v(x) is the realization of a stochastic process. 
In Equation (4), v(x) has the mean zero, variance cr 2 , and non-zero covariance. The weight of 

A 

electromagnet Weight total is replaced by y(x) to make a surrogate approximation model. Let y(x) be an 
approximation model. Hereafter, A means the estimator. When the mean squared error between y(x) 

A A 

and y(x) is minimized, y(x) becomes: 

Xx) = p+r r (x)R- 1 (y-Pq) (5) 

where r is the correlation vector, R is the correlation matrix, y is the observed data and q is the unit 
vector. The definitions of R and r are well explained in Refs [8,9]. 

The unknown correlation parameters of #7, 62,..., 6 n defined in R are calculated from the 
formulation according to: 



maximize - 



[n s ln(G 2 ) + ln\R\] ^ 



where 6 t (i = l,2,...,n) > 0. In this study, a GRG (generalized reduced gradient) algorithm built in the 
Excel program was utilized to determine the optimum parameters. 

4. Optimum Design 

4.1. Metamodel and Design Variable 
Optimization formulation: 

Minimize: Weight tota i(dv h . . ., dv 7 ) (7) 

Subject to F norma i (dv h ...,dv 7 )> 61 1.689 (8) 

Figure 6 shows the prototype electromagnet of the reference model which is analyzed. The 
comparison of the static force obtained from FEM simulation and experimental test is shown in Figure 
7. From the results, the use of FEM is validated, as it can be observed, since the comparison of normal 
force by simulation and experiment test shown in Figure 7 displays good agreement. Selection of the 
design variables is very important setup in optimization procedure. 
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Figure 6. Prototype of the electromagnet. 




Figure 7. Comparison of FEM and experiment (reference model). 
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Seven dimensions are selected as the design variables as shown in Figure 8, which shows the flux 
density vector in electromagnet of 2D and 3D FEM. The magnetization characteristics of the S20C, 
SM490A material is as shown in Figure 9. Table 1 shows the design variables and levels. 



Table 1. Design variable and level. 
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Figure 8. Design variables and flux pattern (2D and 3D FEM). 
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Figure 9. Magnetization curve of the core material used. 
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4.2. Response Surface Methodology 

Table 2 shows the table of mixed orthogonal array and simulation results by 2D FEM. The mixed 
orthogonal array table is determined by considering the number of the design variables and each level 
of them. After getting experimental data by FEM, the function to draw a response surface is extracted. 
In order to determine equations of the response surface for each response value (weight, normal force), 
several experimental designs have been developed to establish the approximate equation using the 
smallest number of experiments. The purpose of this paper is to minimize the objective function 
{Weight tota i) with constraints of normal force {F normal)- The two fitted second order polynomial of 
objective functions for the seven design variables are as follows: 

F normal = 148.75 + 26.1^ - 2.4 ldv 2 -10.02rfv 3 -32.18Jv 4 -l.63dv 5 + 35.16</v 6 

+ 3.0k/v 7 -l.\AE-2dv 2 + \.3E-3dv 2 2 +0.23dv 3 2 + 0.35dv 4 2 + 5.56E-2dv 5 2 
-0.9dv 6 2 -3.12E-4dv 2 

(9) 

Weight total = 8.765 + 0.212^ -0.248dv 2 -0.154</v 3 -0.488dv 4 -0.113rfv 5 + 0.264</v 6 
+ 0.065</v 7 + \ME-3dv 2 + 2.%E-3dv 2 + 6.5AE-3dv 2 + 5.95E-3dv 4 2 
+ 7.58E-3dv 5 2 -4.05E-3dv 6 2 + 2.64E-6dv 7 2 

2 

The adjusted coefficients of the multiple determination R a dj for normal force and weight are 
weighttotal (99%) and F norm ai (100%). Normal force and weight of experiment result according to 
change of the design variables are shown in Table 2. As many parameters are defined as design 
variables, a large simulation time is required due to the large number of required experiments. 
Therefore, it is necessary to establish the significant parameters to investigate the influence on the 
design result. The Pareto chart of normal force and weight shows the magnitude and importance of an 
effect. This chart displays the absolute value of the effects. The geometries of electromagnet can be 
defined by seven parameters, as shown in Table 1 and Figure 8. The most effective design variables of 
normal force and weight are dv 7 , dvj according to the Pareto chart in Figure 10. 




J. j I * 1 i I j L 
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Table 2. Table of mixed orthogonal array Ljs(2 x 3 ). 
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Figure 10. Pareto chart of the standardized effects (alpha=0.05) (a) Normal force response 
(b) Total weight response. 
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The reference line corresponds to alpha = 0.05; 95% confidence interval. 

The values of ineffective design variables are determined by RSM. Figure 11 shows response 
optimization to find optimal solution according to response curves. The slope of the response function 
in Figure 1 1 shows the sensitivity of design variable. According to the response optimization the most 
sensitive design variables of normal force and weight are dv7, dvj. 



Sensors 2010, 10 



6727 



Figure 11. Response optimization. 
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4.3. Kriging Interpolation Method 

The purpose of this paper is to minimize the objective function (Weight tota i) with the constraint of 
normal force (F norma i). The kriging interpolation method, the optimum parameters of 67 are 

determined by solving Equation (6). Then, the estimator P is calculated. In Tables 3 and 4, the optimal 
point is searched to find the point of less than 1 1.412% of the weight and greater than 7.754% of the 
normal force of the initially designed electromagnet in the magnetic levitation system. The optimum 
estimators are shown as Table 4. The optimum solution is the same in the RSM and kriging 
interpolation method. The reason is that the optimum solution is near the boundary value. However, 
response values, Weight to tai and F norma i by RSM and kriging method, respectively, are different. The 
simulation result of the predicted optimum set is shown Table 5, displaying good agreement. 

Table 3. Optimum level and size. 



Design variable 



Model 


dvi 


dv 2 


dv 3 


dV4 


dv 5 


dv 6 


dV7 


Initial 


20 


50 


20 


40 


15 


20 


180 


Optimum (RSM) 


23.877 


45 


16 


40 


7 


16 


166.844 


Optimum (Kriging) 


23.877 


45 


16 


40 


7 


16 


166.844 



Table 4. Optimum P and correlation parameters for kriging models (N s =100). 



Response 




Correlation parameter (corresponding design variable) 




di (dn) 


e 2 {dv 2 ) d s (dv 3 ) 0 4 {dv 4 ) d 5 (dv 5 ) 0 6 (dv 6 ) 0 7 (dv 7 ) 


P 


Weight tota i 


3.026e-3 


2.304e-4 0.909e-3 0.201e-3 6.318e-5 1.195e-3 1.289e-2 


16.0023 


-^normal 


1.319e-2 


1.382e-5 3.437e-6 1.172e-5 7.819e-6 2.320e-6 4.087e-2 


587.009 
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Table 5. Comparison of initial and optimum model. 



Model 



Weight (kg) Normal force (N) 



Initial 



2D FEM 



13.319 578.5 

13.319 573.48 

0 0.86 

12.300 637 

12.107 611.68 

11.799 611.70 

-4.073 -3.972 

-2.610 0.003 

-11.412 7.754 



3D FEM 

Error(2D vs. 3D) % 



Optimum RSM (predicted) 



Kriging method (predicted) 
FEM (verification) 
Error (RSM vs. FEM) % 
Error (Kriging vs. FEM) % 



Variation between initial and optimum FEM % 



5. Conclusions 

This paper deals with optimum design of a lightweight levitation electromagnet on a vehicle, which 
also provides passive guide force, in a magnetic levitation system for contactless delivery applications. 
The optimum design procedure is introduced to design of electromagnet in the magnetic levitation 
system to reduce its weight and to improve the normal force of the initially designed electromagnet in 
the magnetic levitation system using several design variables. The most effective design variables are 
extracted by Pareto chart. The most desired set is determined by RSM and the kriging interpolation 
method and the influence of each design variables on the objective function can be obtained. This can 
efficiently increase the precision of the optimization and reduce the number of experiments in the 
optimization design using the proposed methodologies. 
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